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We propose the Gauss-Bonnet dark energy model inspired by string/M-theory where standard 
gravity with scalar contains additional scalar-dependent coupling with Gauss-Bonnet invariant. It 
is demonstrated that effective phantom (or quintessence) phase of late universe may occur in the 
presence of such term when the scalar is phantom or for non-zero potential (for canonical scalar). 
However, with the increase of the curvature the GB term may become dominant so that phantom 
phase is transient and w = — 1 barrier may be passed. Hence, the current acceleration of the universe 
may be caused by mixture of scalar phantom and (or) potential/stringy effects. It is remarkable 
that scalar-Gauss-Bonnet coupling acts against the Big Rip occurence in phantom cosmology. 



I. INTRODUCTION 

It became clear recently that late-time dynamics of the current accelerated universe is governed by the mysterious 
dark energy. The interpretation of the astrophysical observations indicates that such dark energy fluid (if it is fluid!) 
is characterized by the negative pressure and its equation of state parameter w lies very close to —1 (most probably 
below of it). Quite possible that it may be oscillating around —1. It is extremely difficult to present the completely 
satisfactory theory of the dark energy (also due to lack of all requiered astrophysical data) , especially in the case of 
(oscillating) w less than — 1. (For instance, thermodynamics is quite strange there with possible negative entropy 0). 

The successful dark energy theory may be searched in string/M-theory. Indeed, it is quite possible that some un- 
usual gravity-matter couplings predicted by the fundamental theory may become important at current, low-curvature 
universe (being not essential in intermediate epoch from strong to low curvature). For instance, in the study of string- 
induced gravity near to initial singularity the role of Gauss-Bonnet (GB) coupling with scalar was quite important 
for ocurrence of non-singular cosmology 0, (for account of dilaton and higher order corrections near to initial 
singularity, see also Q). The present paper is devoted to the study of the role of GB coupling with the scalar field 
to the late-time universe. It is explicitly demonstrated that such term itself can not induce the effective phantom 
late-time universe if the scalar is canonical in the absence of potential term. It may produce the effective quintessence 
(or phantom) era, explaining the current acceleration only when the scalar is phantom or when the scalar is canonical 
with non-zero potential. It is interesting that it may also have the important impact to the Big Rip singularity 
similarly to quantum effects 0,0> preventing it in the standard phantom cosmology. Note that we concentrate mainly 
on the exponential scalar-GB coupling and exponential scalar potential, while the consideration of other types of such 
functions and their role in late time cosmology will be considered elsewhere. 

II. THE ACCELERATED UNIVERSE FROM SCALAR-GB GRAVITY 

We consider a model of the scalar field cj> coupled with gravity. As a stringy correction, the term proportional to 
the GB invariant G is added: 

G = R 2 ~ 4R MV RT + R^R^f" . (1) 
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The starting action is given by 

S = J d^x^-g i^R- - V(<P) + /(0)G| . (2) 

Here 7 = ±1. For the canonical scalar, 7=1 but at least when GB term is not included, the scalar behaves as 
phantom only when 7 = — 1 Q showing in this case the properties similar to quantum field 9]. In analogy with 
model [Tof where also non-trivial coupling of scalar Lagrangian with some power of curvature was considered, one 
may expect that such GB coupling term may be relevant for the explanation of dark energy dominance. 
By the variation over <f>, we obtain 

O = 7 V 2 0-y'(0) +/'(</>)£. (3) 
On the other hand, the variation over the metric g pv gives 

= -L (-BT + \gTR) + 7 (~£> M W - \g^d p W^ + \g^ (-V(<t>) + f{<}>)G) 

-2f(^)RR^ + 2VV" (f(ct>)R) - 2 5 ^V 2 (f(cf>)R) 
+%f{4>W p R l/p - 4V P V" (ft<t>)R vp ) ~ 4V P V (/(^iff") 

+4v 2 u(mn + wv p v a {f{m pa ) - v{4>w paT R v Pa r + 4v p v CT {f(m w n ■ (4) 

By using the identities obtained from the Bianchi identity 



V p V a R' ipv ' T = V 2 i?^ - -VVJ! + R ppva R pa - R p p R vp , 
V p V^i? pl/ + V p V v R p i* = - {\7 P \7 V R + V V \7 P R) - 2R ppva R pa + 2R" p R vp 



V p V«R pa - infl, (5) 



one can rewrite Q as 



= ^ + \g» v R) + 7 - ^"dptdPctj + \g p » {-V{ct>) + f(cj>)G) 

-2/(0)iLH^ + ifi^R^P - 2f{$)R^R v p!TT + 4f(<f>)IV>P° v R pa 

+2 (VT/W) i? - 2<T (V 2 /(0)) i? - 4 (V P V7(0)) ii"" - 4 (VpV/M) i? w 

+4 (V 2 /(^)) + 4g p » (VpVa/W) ~ 4 (V„V ff /(<£)) . (6) 

The above expression is valid in arbitrary spacetime dimensions. In four dimensions, the terms proportional to /(</>) 
without derivatives, are cancelled with each other and vanish since the GB invariant is a total derivative in four 
dimensions. 

The starting Friedmann-Robertson- Walker (FRW) universe metric is: 

3 

ds 2 = -dt 2 +a{t) 2 Y,(dx t ) 2 , (7) 



where 



= a 2 Hh 3 , T) t = T l tJ = HS l 3 , Ri tjt = - (h + H 2 ) 6 tJ , R l]kl = a 4 H 2 (fofy - S a S kj ) , 
R tt = -%(h + H 2 ) , Rij =a 2 (h + 3H 2 ) Sij , R = 6H+ 12H 2 , other components = , (8) 



(here the Hubble rate H is defined by H — a/ a). Assuming <f> only depends on time, the (/x, v) = (t, t)-component in 
Q has the following simple form: 

= -^H 2 + U 2 + V^)-24 t pf'(cf > )H 3 . (9) 
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On the other hand, Eq.J2J becomes: 

0= -7(0 + 3i?0) - V'(0) + 24/'(0) (ffff 2 + H 4 ) . (10) 
We now consider the case that V((f>) and /(</>) are given as exponents with the constant parameters Vo, /o, and <f>o 

V = Voe~% , /(</>) = f e% . (11) 

Assume that the scale factor behaves as a = a t h ° ( power law). In case that h is negative, this scale factor does not 
correspond to expanding universe but it corresponds to shrinking one. If one changes the direction of time as t — > — f, 
the expanding universe whose scale factor is given by a = ao(—t) h ° emerges. In this expression, however, since ho is 
not always an integer, t should be negative so that the scale factor should be real. To avoid the apparent difficulty, 
we may further shift the origin of the time as t — > — t —* t s — t. Then the time t can be positive as long as t < t s . 
Hence, we can propose 

H=^, = o lnf, (12) 



when ho > or 



H =--!*- ^ = ^ ln^, (13) 



when ho < 0, with an undetermined constant t\. By the assumption (|12|l or ljl3|) . one obtains 



from @ and 



= 7(1- 3&„) $ + 2Vot\ + (h - 1) , (15) 



H 



from l|T0)l. Using O and (JTHJ, it follows 



48/o^ _ 6 ^ o _7^!V (16) 



*? K 2 (1 + ho) 

The second Ea. fT^|l shows that if — 1 < h < and 7 = 1, / should be negative. Without the GB term, that is, 
fa = 0, a well known result follows: 

h = ^. (17) 

Since the equation of state parameter w is given by 

^ = -1 + ^-, (18) 
3ft 

if ho < (ho > 0), w < —1 (w > —1). Eas. nitty indicate that even if 7 = 1, with the proper choice of parameters ho 
can be negative or w < — 1. Even if 7 > 0, when ho < — 1, Vo is positive, which means that the potential V(0) is 
bounded below. As a special case we consider 



2 = 6h% (1 - fag) 

7(1 -5M« 2 ' 



(19) 



which gives V(cf>) = 0. In order that <fro could be real, one has 



- < ho < 1 , when 7 = 1, or ho > \ or h > 1 . (20) 
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FIG. 1: The qualitative behavior of versus ho from l|19fL 



In the case <|19f) . the scalar field is canonical (7 = 1), and there is no potential V(4>) = 0, even if we include the 
term proportional to the GB invariant, we cannot obtain the effective phantom cosmological solution with ho < or 
w < — 1. Ea. H16(l tells, however, when 7=1 and Vo > even if Vq is arbitrary small, if we choose /o properly, we may 
obtain the effective phantom. The qualitative behavior of 70q versus ho when Vq = is given in Figure^ There is one 
positive solution, which may mimic the effective matter with 1/5 < ho < 1 when 7=1. We also find, when 7 = — 1, 
there are always three solutions for ho from (|19fl . one is given by ho < and describes the phantom cosmology, one is 
ho > 1 describing the quintessence cosmology, and another corresponds to the matter with < ho < 1/5. Then even 
if 7 = — 1, there appear the solutions describing non-phantom cosmology coresponding the quintessence or matter. 
As an example, we consider the case that 



80 



-1 



which gives, from (|18|) . 



-1.025 



(21) 



(22) 



This is consistent with the observational bounds for effective w (for recent discussion and complete list of refs., see[rH|). 
Then from l|16|) . one obtains 



2 1 ^ 531200 



403 
154 



70o «' 



Jo 



1 / 9 



27 



k 2 V 49280 7884800 



70o 



(23) 



Therefore even starting from the canonical scalar theory with positive potential before introducing the term propor- 
tional to the GB invariant, we may obtain a solution which reproduces the observed value of w as in Ij22(l . 
In case of the model induced from the string theory 0, we have Vq = (V(0) = 0) and 



in lfTT|) . Then Eq. ljT§|l reduces as 



which has only one real solution as 



3hl - 3hi + 5h - 1 = , 



h = 0.223223 



The solution gives 

w = 1.98654 . 

There is another solution of © and l|10|) with i|llf) . In the solution, and H are constants, 



(24) 



(25) 



(26) 



(27) 



= <A) , H = Ho , 



(28) 
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what corresponds to deSitter space. Using © and l|10[l with (|llfl . one finds 

_ 2< fo 

e *o 

ff o--^^- (29) 

Therefore in order for the solution to exist, we may require fo < 0. In (|29|l . <po can be arbitrary. Hence, the Hubble 
rate H = Hq might be determined by an initial condition. 

In case of the model , the term including the GB invariant always gives the contribution in the same order with 
those from other terms even if the curvature is small. This is due to the factor /(</>), which enhances the contribution 
when the curvature is small. 



III. LATE-TIME ASYMPTOTIC COSMOLOGY IN SCALAR-GB GRAVITY AND BIG RIP 

AVOIDANCE 

In the following, another model, which is slightly different from may be considered: 

V{()>) = V e~% , /(</>) = / e^ , (a > 1) , (30) 

Different from the model , the model l|3UI) will not be solved exactly. We can only find the asymptotic qualitative 
behavior of the solutions. Nevertheless, the asymptotic behavior suggests the existence of the cosmological solution, 
where the value of w could vary with time (oscillation) and/or could depend on the curvature. 

Assuming the solution behaves as i|12|) or (|13|l . when the curvature is small, that is t in (|12|l or t s — t in (|13|l is 
large, the GB term becomes small and could be neglected since it behaves like l/t~a+ 4 or 1/ (t s — t) ~ = +4 . When the 
curvature is small, the solution could be given by (|17|l . then the effective phantom phase with w < — 1 could appear 
only in case 7 = — 1 < 0. On the other hand, when the curvature is large, that is t in (|12|) or t s — t in (|13|l is small, 
the classical potential could be neglected. Without the classical potential, by assuming, instead of (|12|l 

H = ^ , <j> = acpa In , (31) 



when ho > or 



t s -t 

when ho < 0, the following equations replace (|14|l and (|15(l 



H = < f, = o^ a hi^—^ , (32) 







.2 



= 7 (l-3/ l0 )« 2 ^ + ^ 2 ^(/' ! » ■ 1) ■ (•'«) 
By deleting fo in the above two equations, one gets 



6/ig (1 - h ) 



ja 2 (1 — 5ho) k 2 



(34) 



which corresponds to (|19|l . Then when 7 = 1, the solutions of Q34f) are not qualitatively changed from those of i|19fl . 
and there is only one solution 1/5 < ho < 1. On the other hand, when 7 = —1, since the sign of the r.h.s. in l|19(l 
is changed from 7 = 1 case, as clear from FIG^ there are three solutions, corresponding to the phantom ho < or 
w < —1, the quintessence ho > 1 or —1 < w < —1/3, and the matter with < ho < 1/5 or w > 7/3. Then if the term 
proportional to the GB invariant in case 7 < (which corresponds to a scalar phantom solution without GB term) is 
included the effective w can become larger than —1 and the Big Rip singularity might be avoided (see p] for quantum 
effects account to escape of Big Rip). That is, in case 7 < 0, when the curvature is small as in the current universe, 
the GB term is negligible and the potential term dominates, which gives the cosmic acceleration with w < — 1. Then 
the curvature increases gradually and the universe seems to tend to the Big Rip singularity [{|- However, when the 
curvature is large, the GB term becomes dominant and might prevent the singularity. Hence, in case 7 < 0, the GB 
term may work against the Big Rip singularity occurence, like quantum effects Q. After the GB term dominates 
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when 7 < 0, the curvature turns to become smaller. Then the potential term dominates again. This might tell that 
the behavior of the universe might approach to the deSitter sapce with w — — 1 by the dampted oscillation. In fact 
even in the model (pTTTft - if (|2"%)> is assumed, there is a deSitter solution corresponding to 



In (|3Tj|) . we have assumed a > 1. If we consider the case that V{4>) — Vbe~*o and /(</>) = / e°*« as in (|3"Ujl but 
< a < 1, there appears a solution where the term including GB invariant becomes dominant even if the curvature 
is small. By assuming l|31|) or Ij32|l . Eq.(J2HJ) is obtained again. Hence, when < a < 1, the solution where ho can be 
positive or w > -1 even if 7 < (scalar phantom) appears. 

On the other hand, if 7 > there is no accelerated universe solution with w < — 1. The parameter w may change 
with time but w is larger than —1/3. curvature due to the GB term, effective equation of occur eventually. 



We considered essentially two models with exponential couplings given by (| 1 1 ft and (|30|l . The model (|11J) may be 
considered as the special case corresponding to a = 1. The main results can be summarized as follows: 

1. a = 1 case: exactly solvable 

(a) V = case: When 7 = 1, there is only one solution —1/3 < w < 7/3. On the other hand, when 7 = — 1, 
there are three solutions, corresponding to w < —1, —1 < w < —1/3, and w > 7/3. 

2. a > 1 case: the potential term dominates for small curvature and the GB term for large one. 

(a) 7 > 0: The value of w may be time dependent but there is no solution describing acceleration of the 



(b) 7 < 0: There might appear the Big Rip singularity but there might be a solution asymptotically approaching 
to the deSitter space. 

3. < a < 1 case: the potential term dominates for large curvature and the GB term for small one. 

(a) 7 > 0: There is no solution describing acceleration of the universe. 

(b) 7 < 0: There appears the Big Rips singularity. 

For the models (|11[) and H3U|) . in case Vo = (that is, when the potential vanishes), by replacing acf>o with 0o, it 
follows the two models are equivalent. Especially 7 = — 1 case has been well studied and it has been shown that there 
are always three effective cosmological phases corresponding to the phantom with ho < or w < — 1, the quintessence 
with ho > 1 or — 1 < w < —1/3, and the matter with < ho < 1/5 or w > 7/3. Even if Vq 7^ 0, the model i|ll|l can be 
solved exactly and the solutions where ho and therefore w are constants may be found. On the other hand, if Vo =/= 0, 
in the model (|30() there exist the solutions where the values of ho and therefore of w are time-dependent. There could 
emerge a cosmology, which behaves as phantom one with w < — 1 when the curvature is small and as a usual matter 
dominated universe with w > — 1 when the curvature is large. Moreover, Big Rip singularity does not occur. 

Our study indicates that current acceleration may be significally influenced by stringy/M-theory effects (terms) 
which somehow became relevant quite recently (in cosmological sense). It remains a challenge to construct the 
consistent dark energy universe model from string/M-theory. 

Note added: After the first version of this paper (with some error) appeared in hep-th, the related study, for 
instance, of the influence of scalar-GB term to Big Rip has appeared in ref . jl2| . 
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IV. DISCUSSION 



universe. 
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APPENDIX A: STABILITY OF PHANTOM COSMOLOGY 

In this appendix, we check the stability of the above solutions. The following quantities are convenient to introduce: 

„,, ,. d d 1 d , . . 

For simplicity, we also put k 2 to be unity. Then by using ijjjj and i|l(J|) with one finds 



dX 



7 2 O X 3 + 2-fX (8X 2 Z - O (3 + 52XZ)) + 4 (Zgjts. + 24V ° f ° X + 12Z (0 O + 160 O XZ - 8X 2 Z)^ 



N 20 o (7 + 6jXZ + 96Z 2 ) 

dZ Z (- 7 2 0o* 2 - 1GZ + 12 (0 O - X) Z) + 2 7 (X + 160 O XZ)) 

dN ~ O (7 + 6jXZ + 96Z 2 ) 

For the solution (JTIJJl or JTSJl , it follows 



(A2) 
(A3) 



X = X ^^ , Z = Z ^^. (A4) 



In terms of Xo and Zo, Eas. (|14|l and (|15|) can be rewritten as 

= -4 + ^ + ^-24^0, (A5) 
= 7 X 2 - 30o7*o + + 240o^o - 2AZ X Q . (A6) 

00 A) 

For the solution l|A4|l . by using (| A5|) and (|A6|I . the right hand sides of Eas. l|A2|) and (|A3|) vanish consistently. We 
now consider the perturbation around the solution l|A4fl : 

X = X a + SX , Y = Y + SY . (A7) 

We now only check the stability for V = (Vb = 0) case. 
Using (|A2|) and (|A3|) . one obtains 



Here 



.4 
B 
C 



37 2 0qX o 2 + 487X5% - 6 7 0q - 2O8 7 0oX o Z o + 768 (0 O - X ) Z 2 

20o ( 7 + 8jX Z + 96Z 2 ) 
16 7 X3 - 1O4 7 O X 2 + 480o + 15360 O X O Z O - 768X6% 

20o (7 + 8-fXoZo + 96Z 2 ) ' 
2Z (-7 2 0qX o + 96Z 2 + 7 + 16 7 0oZo) 
0o (7 + 8 7 X Z + 96Z 2 ) 



b 32Zo(-12(0o-Xo)Zq + 70oXo) 

0o (7 + 87X0Z0 + 96Z 2 ) " 1 ' 

If the real parts of all the eigenvalues of the matrix M are negative, the perturbation becomes small and the system 
is stable. Then the condition of the stability is given by 

A + D<0, AD-BC>0. (A10) 

First, the check of the stability for V — case is in order. By using (|A4|) . JUJl, and (|A6|) . we find 



2 _ 6 (h - 1) 2 _ 7 (3feo - 1) 3/ip-l 

h 2 7(5^0-1)' ° 96 (/10 - 1) (5/10 - 1) ' °° 4(5/io-l) 
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In order that Xq and Zq are positive, it follows 

| < h < 1 , when 7 > 

or ho < \ or ho > 1 , when 7 < . (A12) 

By using l|Allj) . A, B, C, and Z? in (|A9|) can be expressed in terms of ho: 

- _ (fop - 1) (9/ig - ihp + 1) - _ 24 (ftp - 1) (3fe 2 - 2h Q + 1) 

ho (5/ig - 4/i + 1) ' - 7/io(5^-4ft + l) 

- 7 (3fto-l)(3feg-l) - 2 (h - l) 2 (3fe - 1) 

12 (/i - 1) (5fc§ - 4fc§ + i) ' ft,, ( 5 ft2 _ 4/10 + 1) ' 1 j 

Then we obtain very simple results: 

A + D = - 3 J^1, (A14) 

AD-M= 2(1 -f o) . (A15) 
h 

Therefore Eq. (|A10fl is satisfied and the system is stable if and only if 

ho < . (A16) 
Then the case corresponding to phantom cosmology with ho < is always stable. 
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